Weak localization and magnetoconductance of Dirac fermions under charged 

impurities in graphene 
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On the basis of self-consistent Born approximation, we present a theory of weak localization of 
Dirac fermions under finite-range scatters in graphene. With an explicit solution to the ground state 
of singlet pseudospin Cooperons, we solve the Bethe-Salpeter matrix equation for all the singlet and 
triplet pseudospin Cooperons at long-wave length states by perturbation treatment. The solution to 
the Cooperon in the presence of the external weak magnetic field is also obtained. We calculate the 
quantum interference correction to the conductivity and present the comparison with experiments. 
It is shown that the present calculation for the magnetoconductivity is in good agreement with some 
of the experimental measurements. 
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I. INTRODUCTION 



It has been found that the charged impurities with 
screened Coulomb potentialsi*^!^ are responsible for the 
observed carrier density dependence of the electric con- 
ductivity of graphene.— In a recent work, we have investi- 
gated the weak localization (WL) of electrons under the 
charged impurity scattering in graphene.^ The descrip- 
tion for the Cooperons under the finite-range scatters is 
different from that for the zero-range potentials as stud- 
ied in the existing works,-^'^>S In this paper, we present 
the details of the formalism for the WL of Dirac fermions 
under finite-range scatters in graphene. We also calculate 
the quantum inetrference correction (QIC) to the electric 
conductivity under a weak magnetic field and compare 
the result for the magnetoconductivity with the experi- 
mental measurements. 

The central problem of theoretically studying the weak 
localization of Dirac fermions under finite-range scatters 
is to solve the Bethe-Salpeter matrix equation for the 
Cooperons. With the self-consistent Born approximation 
(SCBA) to the single particle, we can obtain an explicit 
solution to the ground state of the singlet pseudospin 
Cooperons. By perturbation method, we will solve the 
Bethe-Salpeter matrix equation for all the singlet and 
triplet pseudospin Cooperons at long- wave length states. 
With the Cooperons, we derive the quantum interference 
correction to the electric conductivity that gives rise to 
WL effect. 

At low carrier doping, the low energy excitations of 
electrons in graphene can be viewed as massless Dirac 
fermions ^liiii^iii^'i^ This has been confirmed by recent 
experimentsi^ii^ Using the Pauli matrices cr's and t's to 
coordinate the electrons in the two sublattices (a and b) 
of the honeycomb lattice and two valleys (1 and 2) in 
the first Brillouin zone, respectively, and suppressing the 
spin indices for briefness, the Hamiltonian of the system 



is given by 



k kq 



(1) 



where "01 



(4 



feal' '-febl' ^kb2^ 



4a2) t'^^ fermion operator, 
the momentum k is measured from the center of each val- 
ley, V (~ 5.856 eVA) is the velocity of electrons, V is the 
volume of system, and Vi (q) is the finite-range impurity 
potential. For charged scatters, Vi{q) is given by. 



n'i.{-q)vo{q)cro n,{Q - q)viai 
nii~Q ~ q)viai ni{-q)vo{q)ao 



(2) 



where ni{~q) is the Fourier component of the impurity 
density, VQ{q) and vi are respectively the intravalley and 
intervalley impurity scattering potentials, and Q is a vec- 
tor from the center of valley 2 to that of the valley 1 [Fig. 
1(a)]. In Appendix, we detail the discussion on this im- 
purity potential. Here, all the momenta are understood 
as vectors. 

Under the SCBA [Fig. l(b)]^iiiii^ the Green function 
G{k,uj) = [uj-^- ^ — vk-aTz~Ti(k,uj)]~'^ and the self-energy 
Yi{k,uj) = YiQ{k,uj)'\-Y,c{k,ijj)k-aTz of the single particles 
are determined by coupled integral equations)^ 



So(fc, w) 

Sc(fc,Cj) 



V 
V 



Y,[vl{\k-k'\) + vl]G,{k',u) (3) 
fe' 

Y,vl{\k-k'\)G,{k',w)k-k' (4) 



with 



hk 



Go{k,u) 
Gcik,uj) 



where uj = + — So(fc, w) with /i the chemical potential, 
hk = vk+T^c^k, uj), k is the unit vector in k direction, and 
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FIG. 1: (color online) (a) Brillouin zone and the two Dirac- 
cone valleys, (b) Self-consistent Born approximation for the 
self-energy. The solid line with arrow is the Green function. 
The dashed line is the effective impurity potential, (c) Cur- 
rent vertex with impurity insertions. 

the freqtiency lo is understood as a complex quantity with 
infinitesimal small imaginary part. The current vertex 
vTx{k,ijJi,uj2) [Fig. 1(c)] can be expanded as 

3 

r^(fc,wi,W2) = ^j/j(fc,wi,W2)A;^(fc) (5) 
where Ag(fc) — t^cTx, Af{k) = (Jx(^ ■ k, Af(fc) = (? ■ 



kax, A^ik) = TzG ■ kax(J ■ k, and yj(k,u!i,uj2) are de- 
termined by four-coupled integral equations;^ The x- 
direction current-current correlation function [Fig. 2(a)] 
is obtained as 



P(a;i, W2) = 




with X,{k,L0i,L02) = Tr[G(fc,c^i)^J(fc)G(fc,a;2)^g(fc)], 
for w's (wi and L02) = w ± iO = w^. The detailed deriva- 
tions of ra;(A:, wi, 0^2) and P(ijji,uj2) can be found in Ref. 
and will not be repeated here. 



II. FORMALISM 

The WL effect in the electric conductivity stems from 
the QIC to the electric conductivity. Theoretically, it is 
given by the maximum crossing diagrams as shown by 
Fig. 2(b)iiii^ii2ii^ The process of the maximum cross- 
ing diagrams is associated with two-particle propagator 
Ca^aia^aS^.k',q,Lo) (Coopcrou). It obcys the Bethe- 
Salpeter 16 x 16 matrix equation represented in Fig. 2(c). 
Here, the superscripts j's denote the valley indices, and 
the subscripts a's correspond to the sublattice indices. 
To explicitly write out the equation of Fig. 2(c), we here 
give the simpler one for Ctiaiaaoni^^ ^'j 1^ ^) that starts 
from the single impurity line [the dashed line with a cross 
in Fig. 2(c)], using the convention J^'^^-'^^ — Sj^j^Sa^ai and 
j (a) as the conjugate valley (site) of j (a): 



+ y E "o(|fc - ^il)GlV,^(9/2 + fci,c.+)Gi^/J,(9/2 - fci,c.-)C^;3',=i«^^(fc^ k',q,u;) 

fci/3/3' 

+ y E ^lGY^0(<l/'2 + k^,u;+)G^^^^,,{q^ (6) 



The first term in the first line in right hand side of Eq. 
([6|) is due to the intravalley scatterings, while the second 
term comes from the intervalley scatterings. Sj_^j^ means 
that when a particle is scattered to valley ji, another 
particle should be scattered to valley ji so that the to- 
tal momentum (vanishing small under consideration) of 
the Cooperon is unchanged. The second and third lines 
are the processes of Cooperon propagating after the in- 
travalley and intervalley scatterings, respectively. The 
equation for C is obtained from Eq. ^ by subtracting 
the single impurity line from C. 



The form of Eq. ([6]) seems rather miscellaneous. It 
may be simplified by classifying it with good quantum 
number of the Cooperons. To do this, we note that the 
elements of the coefficient matrix of C in Eq. ^ are 
arranged according to the indices (superscripts and sub- 
scripts) of the Green functions. Since the Green function 
G(fc, Lu) are composed by the unit matrix and r^t? • k, we 
then look for all the operators that commute with t^ct- k. 
Recently, McCann et al.^ have introduced the operators 
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FIG. 2: (color online) (a) Electric conductivity, (b) Quan- 
tum interference correction to the conductivity, (c) Cooperon 
propagator. 

of isospin S's and pseudospin A's, 

So =To(To, ^x=Tzax, T,y=Tz(7y, Y.z~Toaz, 
Ao —ToCFo, Kx—Tx<Jz, ^y—TyUz, Az=Tzao- 

Clearly, A's commute with E's and TzCT • k, and are con- 
serving operations for the Cooperons. Therefore, we 
transform the Cooperons from the valley-sublattice space 
into the isospin-pseudospin space according to McCann 
et ai, 

pll' _ }_ fM'yi^^ fhhhu (7\ 

where Mj = SyE^AyA;. We will hereafter occasionally 
use the indices 0,1,2,3 or 0,x,y,z to label I and s. In the 
isospin-pseudospin space, the single impurity line is given 
by Wi'^, = W^Si^l, with 

WU\k - fell) = n,vl{\k - fell) + n,vl{Sio - <5^.)(-l)^ 

which is the transform of the first line in the right hand 
side of Eq. ([6]). The result of second-|-third lines in right 
hand side of Eq. (5]) is transformed to 

^ [W'i\k - k^\)hiki,q)UCi'lAkl,k', q) 

where h{ki,q) is a matrix defined in the isospin space 
with the element given by 

hss'{ki,q)=TiiG{-k+,Lj+)J:sG{^k^,Lo-)j:l,]/4 (8) 

and kf = ki ± q/2. With these results, we obtain the 
equation for Cf^, , 

C'jUk,k',q) = 1 ^nL^(fc,fci,g)[W^i^(|fci-F|)C' 
+ Cll,ikuk',q)] (9) 



where Il''{k,ki,q) = W'{\k — ki\)h{ki,q). Here, the ar- 
gument LU of C"^, and II' has been suppressed for brief- 
ness. From Eq. ([H]), it is seen that the pseudospin of the 
Cooperon is indeed conserved during the impurity scat- 
terings, C"g, = C^g, Su> . We then need to deal with C^^,, . 
Thus, the original 16 x 16 matrix equation is separated 
into four 4x4 ones, each of them corresponding to a def- 
inite pseudospin I. In the isospin space, the Cooperon of 
a pseudospin Hs a 4 x 4 matrix denoted as CK 

To solve Eq. we use the standard method that 

expands C' in terms of the eigenfunctions of 
n'(fc,fci,g): 

k', q)^Yl Q^ik', q), (10) 

n 

where cj^(^) are constants and 

i n'(fc, fci, g)*„(fci, g) ^ \Uq)'fn{k, q) (11) 

ki 

with Ajj(q) the eigenvalue. Here, 'i>l^{k,q) is a column 
vector with four components in the isospin space since 
H'(fc, ki,q) is a 4 X 4 matrix in this space. The constants 
c5j(g) are determined by substituting Eq. (fTO|) into Eq. 
([^ . It is then seen thatc'„((7) cx [1-A5j((?)]-i. Therefore, 
the predominant contribution to C' comes from the state 
with the lowest |1 — Ajj(g)| = |1 — A'(q)| that can be 
vanishing small. We will here take into account only the 
state of the lowest \1 — X\q)\ for each /. 

Firstly, we consider the case oi I — and q = 0. A 
solution can be explicitly obtained as A''(0) ~ 1, and 
^'•'(fc,0) = *(fc) with 

**(fc) = [Ao(fc,cj),-Ac(fc,w)cos(/), -Ac(fc,w)sin(/),0] 

where ^*(fc) is the transpose of ^'(fc), Ao(fc,a;) = 
ImEo(fci^ ): Ac(fc,w) = ImEc(fc, w"), and is the an- 
gle of k. The four components of ^'(fc) correspond to 
s = 0,x,y,z respectively. The solution of A"(0) = 1 is 
the most important one which gives rise to the diverging 
contribution to the Cooperon. One may check this result 
with the help of Eqs. ^ and (j!]). Actually, the above so- 
lution is just a consequence of the Ward identity (under 
the SCBA for the self-energy): 

ImE(fc,c^-) = -LYmk-k')lmG{k',Lj-)V,{k' -k)) 

k' 

where (• • • ) means the average over the impurity dis- 
tributions [Fig. 1(b)]. There are three non- vanishing 
components in ^'(fc) because of the finite-range impurity 
scatterings. For the zero-range potential, only the first 
component of \E'(fc) survives and is a constant. One then 
needs to solve a scalar equation instead of the matrix 
integral equation. 

For finite but small q, by expanding lP(k,k',q) to 
second order in q and regarding the difference from 
Il'^(fc, A:', 0) as a small departure, we then solve the 
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FIG. 3: (color online) Eigenvalue X''{q) as function of q at 
electron doping concentration S — 2.0 x 10"'*. 

problem by perturbation method. Since expanding 
II°(fc,fc',(7) [equivalent to expanding h{k',q)] is an ele- 
mentary manipulation but tedious [because there are 16 
elements in h{k',q)], we here just present the result. For 
I ^ 0, the difference between II' and 11° comes from the 
intervalley scattering term in W" . Similarly, we can treat 
this difference by perturbation. For all the cases, to the 
first order in the perturbation, we have 

« X\0)-diq^, for q^O (12) 

where (\E'|5') = J2k '^Hk)'^{k)/V, A'(0) and di are posi- 
tive constants. To the 0th order, the eigenfunction '^{k) 
is unchanged. We here consider only the case of small q 



since that is where QIC is significant. 

In Fig. 3, the eigenvalues A'(g) are shown as functions 
of q at electron doping concentration S = 2.0 x 10~* 
(the doped electrons per site). The impurity scattering 
potential is given by the screened Coulomb one and the 
impurity concentration is chosen as rii = 1.15 x lO^'^a^^ 
(with a the lattice constant).'^ The eigenvalue for I = 
x,y is degenerated. In the limit q — > 0, only 1 — X'^{q) 
approaches zero. The finite value 1 — A'(0) for / ^ is 
determined by the intervalley scattering strength vi. At 
vi = 0, A'(0) = 1 is valid for all LI — A'(0) is larger for 
stronger vi. 

The state for each I so obtained is of the lowest 1 1 — 
A'((7)|. For the lowest state, c^(g) is given by c'((?) — 
q/[1 - \\q)] with 

ci - /^lizya E ^Hm\k,k,,0)W\\k^-k'\)M>{k') 

^ ' ' kkik' 

where the g-dependence of ci has been neglected because 
of the drastic behavior of 1/[1 — A'(g)] at small q. The 
Cooperon is finally approximated as 

C\k,k',q) ^ Ci^{k)^\k')/[1 ~ X'{q)]. (13) 

For the zero-range scatters, only the isospin-singlet 
Coo(fc,fc','?) survives and is independent of k and k' . 
In this case with vi — and lu — 0, hy the one- 
band approximation, one obtains X''{q) = 1 — Dq^ /AAq 
where D — w^/2Ao is the diffusion constant and Ao is 
the first component of 5*. The resultant Cooperon is 
AAoUiVQ/ Dq^ consistent with that of Ref. ^to the order 
of q~'^ in g 0. 

With the Cooperon C', the QIC to the current-current 
correlation function SP{lu^ jLU'^) is calculated according 
to Fig. 2(b). Because the vertex, the Green functions 
and the Cooperons are matrices, one cannot write out 
6P{u!~ ,U!~^) immediately. For doing it, we start to work 
in the valley-sublattice space. According to the Feynman 
rule, we have. 



6P{co~,u;+) = 1^ E[V;(-fc,c.-,^+)]S'/^^C^;i^f^^„^(-^^ (14) 

kqja 



where Vx{k,uji,uj2) = G{k,u;i)rx{k,uii,uj2)G{k,uj2) is 
the vertex connected with two Green functions. With 
the inverse transform of Eq. ([7|) 



4 ^ 

IV 



we obtain 



(5P(w-,w+) 



2V' 



^Tr[Z'(fc,L.)C'(-fc,fc,g)] (15) 



kql 



where Z^{k,uj) is a matrix with elements Z^^, defined as 
Zl, =Tr[V^{k,u;+ ,iu~)M%{-k,u;- ,iu+)Mlr]. 
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There is a simple relation, Z\k, uj) = —Z^{k, us) for I ^ 0, 
because 



!viMlk\ = -Af °, for l = x,y,z 



and the operator A; commutes with G and Tx- This re- 
sult means that the QIC by the pseudospin singlet {I — 0) 
is negative, while it is positive by the pseudospin triplets 
{I = x, y, z). Substituting the results given by Eqs. (fT^ 
and (jl3p into Eq. (|15p and carrying out the (/-integral, 
we get 



ci ,^ 1 - X\0) + diql 



(16) 



r, Ot r. 



(a) 



r, or r. 




'"2 -"^i ''4 



(b) 



FIG. 4: (color online) (a) Cooperon propagator in real space, 
(b) Quantum interference correction to the conductivity. 



where Nq ~ —1, Ni^x,y,z = 1; 9o s-^d 9i the lower 
and upper cutoffs of the g-integral. 

The lower cutoff qg is given by qo = ma,x{L^^^ , L~'^) 
where Lm is the length the electrons diffuse within an 
inelastic collision time Ti„ and L the length scale of 
the system.— At very low doping {S < 1.0 x 10~^, 
the doped electrons per site), Tm due to the inter- 
electronic Coulomb interaction is estimated as Tin ~ 
0A62v/aT'^ (where T is the temperature) from the re- 
cent study of the interacting electrons in graphene us- 
ing renormalized-ring-diagram approximationii^ Li„ is 
then given by = {v^TTm/'^y^'^ where the elastic colli- 
sion time T is determined by the non QIC-corrected con- 
ductivity (Jo, T = hwao/vkFe^ (with hp as the Fermi 
wavenumber).— For low carrier density, we find that Li„ 
is about a few microns for 4 K < T < 20 K. On the other 
hand, the upper limit is qi = Lq ^ with Lq = vt as the 
length of mean free path. 

Using the present formalism, we have recently stud- 
ied the WL effect of Dirac fermions in graphene.^ It 
is found that WL is present in large size samples at 
finite carrier doping. The strength of WL becomes 
weakened/quenched when the sample size < Lin (about 
a few microns at low temperatures) as studied in the 
experimenti^i^i Close to region of zero doping, the system 
may be delocalized. Physically, at small electron doping, 
the Fermi circle and the typical momentum transfer q 
(~ 2kp) are small and also the screening is weak, leading 
to stronger vo{q) than vi. For weak intervalley scatter- 
ing, all A'(0)'s {I ^ 0) close to 1, the QIC from each 
pseudospin channel has almost the same magnitude. Af- 
ter one of ^ 7^ is canceled by the I = channel, the 
net QIC is positive. This is consistent with the fact that 
Dirac fermions cannot be scattered to exactly the back- 
wards direction in case of z;i = and the WL is absent. 
On the other hand, with increasing electron doping, the 
strength of the intervalley scatterings becomes stronger, 
leading to the appearance of WL in large size samples. 
The detailed numerical study of the WL in graphene has 
been presented in Ref. 5 and will not be repeated here. 



We here concisely explain why pseudospin singlet 
Cooperons give rise to WL (negative QIC) but the pseu- 
dospin triplets result in anti-WL (positive QIC) . As seen 
from the matrix Afj defined below Eq. ([7]), the pseu- 
dospins are actually associated with the matrices Aj^A;, 



Aj^Aq = T2(T3, singlet 

Ay Ax ^ —iT3, triplet with I = x 

Ay Ay = 1, triplet with I = y 

AyAz = iTia^, triplet with I = z. 

In terms of the c-operator, 'ip*^^2-k-^y-^i'4'q/2+k annihilates 
a Cooperon of total momentum q and relative momen- 
tum k with pseudospin /. It is seen that in the pseudospin 
singlet state the two particles are in different valleys and 
the parity is odd under the exchange of the valley indices. 
The magnitude of the wave function of the pseudospin 
singlet Cooperon is large when the two particles occupy 
respectively the opposite momentum [defined respect to 
the origin of the Brillouin zone, see Fig. 1(a)] states of 
the single particles. It implies a strong backward scat- 
tering for the electrons and thereby leads to WL. On the 
other hand, the pseudospin triplets are even under the 
valley exchange. For / = z, though the two particles are 
in different valleys, the wave function of the Cooperon 
is small when the two particles occupy respectively the 
opposite momentum states. In this case, the backward 
scattering is weakened, resulting in the increase of the 
conductivity. For / = x and y, the two particles occupy 
the states in the same valley and their total momentum is 
finite. The case corresponds to the final state of the scat- 
tered electrons is not in the backward direction, giving 
rise to a positive contribution to the conductivity. The 
anti-WL can be reduced when the backward scattering is 
strengthened by the intervalley scatterings. 
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III. MAGNETOCONDUCTIVITY 

We here consider that the system is acted with an ex- 
ternal magnetic field B perpendicular to the graphene 
plane. The WL of Dirac fermions in the existence of 
weak magnetic field can be treated in a way parallel to 
Ref. '23. Here, we outline the main steps. Since the sys- 
tem becomes an inhomogeneous one in this case, we need 
to start with the description in real space. The kinetic 
part of the Hamiltonian Hq is 

Ho ^ J drtp''{r)vT^a ■ {p + A)4>{r) (18) 

where p is the momentum operator, and A is the vector 
potential. Here, we have used the units of e = c = 1 
(with — e as an electron charge and c the light velocity). 
With gauge transform 

i^{r) = i^{r) exp(— i / dl ■ A), 

J ro 

I 



A is eliminated from Hq. On viewing this gauge trans- 
form, for very weak magnetic field i?, the Green function 
G{r,r',T — T') — —{Til!{r,T)4'''{r',T')) can be approxi- 
mated as2^ 



G{r,r',oj)^G{r-r',oj)exp{i dl- A), (19) 



where G{r — r',uj) is the Green function in the absence 
of the magnetic field. Here, the position r is understood 
as vector. 

In real space, the Bethe-Salpeter equation for the 
Cooperon is diagrammatically shown in Fig. 4(a). One 
can then write out it explicitly and do the same transform 
as from Eq. ^ to Eq. The final matrix equation 

(in the isospin space) is given by 



C''(ri,r2,r3,r4) = W\ri - r2)[h{ri,r2,r3,ri)W'- {rs - r^) + / dxi / dx2h{ri,r2,Xi,X2)C'^{xi,X2,r3,r^)] (20) 



h°,,{R- R' + ^—^,R- R' - ^—^)exp{i2 dl ■ A) (21) 



where W^{r) is the real space representation of W^{q), and the element of h{ri^ r2, r^, 7*4) is given by 

hss'{ri,r2,r3,r4) = iTr[G(ri, ra, tj+)I]sG(r2, r4, w~)E^,]. 
Using the approximation given by Eq. (|19p . we have 

dl ■ A + i I dl ■ A) 

^ _ J,' r ~ r' 

——,R-R' — 

2 2 

where R — {ri + r2)/2, r = ri — r2, R' = (r^ + 7'4)/2, r' = — r^, and hP is defined in the absence of A. The Fourier 
transform of h^{R - R' + R - R' - ^) is given by 

h\R -R' + ^— ^, i? -R'- ^— ^) ^('^' g)e*« (^-^')+'^ (^"-") (22) 

^ V 

kq 

where h{k,q) is defined by Eq. ([8|). The eigenvalue problem of Eq. (f20ll reads 

W\n~r2) J dr[J dr'2h{n,r2,r[,r'2Mr[,r'2) ^ E',p{n,r2). (23) 

where "0(^11 '"2) is a four-component vector in the isospin space, and is the eigenvalue. Using the coordinates R 
and r, we separate ip{ri,r2) as •4'{ri,r2) — <I>(i?)^(r) where <&(i?) is a scalar representing the motion of center of 
mass, and ^^(r) is a four-component vector meaning the relative motion of the Cooperon. Since the magnetic field is 
weak, only the large-scale motion of the center of mass is significantly affected; the magnetic filed infiuence on the 
relative motion is negligible. Then, '^{r) can be considered as the real space representation of \E'(fc) given in Sec. H. 
By integrating out the relative motion, we get 



dR'>^\q) exp[iq- (i? - R') + i2 J dl- A]<^>{R') = E'<i> 



{R) (24) 



with \\q) defined by Eq. 1^. Using A'(g) « A'(0)-rf;g2 foj. gj^afi q, $(i?') = exp[(^' - R) - V]$(i?), and carrying 
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out the g-summation and i?'-integral, we obtain 

[X'{0)~di{P + 2Af]^R)^E^<i>{R), (25) 

where P = — iV is the momentum operator (of the center 
of mass) of the Cooperon. Using the Landau gauge A = 
(0, Bx, 0), one has 

< = A'(0) -4d,B(n + l/2), (26) 

and $„(i?) being the wavefunction of the corresponding 
Landau state. The degeneracy of each level is g ~ BV/tt. 
The Cooperon is obtained as 

n ^ 

The QIC to the conductivity in the presence of mag- 
netic field is calculated according to Fig. 4(b). To explic- 



itly write out the expression, we note that the magnetic 
field effect on the vertex is negligible small. The strong 
dependence of B in QIC comes from the Cooperon due 
to the small denominator 1 — E'^. Therefore, the vertex 
associated with x' in Fig. 4(b) connected with two Green 
functions after integrating over x' is given by 

i^x' (^3 - ^2 , , c^- ) = E (fci , ^+ , ^^y'"' ' ^"^^ ^ 

ki 

where 14 (fei,w''",w~) has appeared above Eq. ^T5\i . Sim- 
ilarly, for the vertex associated with x in Fig. 4(b), one 
gets 

T4(r4-ri,c^-,c^+) = lEl4(fc2,c^-,c^+)e'"'-(^"^-^~^). 

fc2 

The QIC to the current-current correlation function is 



SP{cj-,u;+) = ^J2 /^''i---/ Mr[Z'(r3-r2,r4-ri,u;)C'(ri,r2,r3,r4)] (27) 

with Zl^,{r3 - r2,r4 - ri,uj) = Tr[y^*,(r3 - r2,uj+ ,uj^)MlVx{ri - ri,LU^ ,lu+)M!.T]. Substituting the results for Vx' , 
Vx, and C'(ri, r2, r3, r4) into Eq. ((27l) . setting ki.2 — ±fc + q/2 and neglecting the small g-dependence in Vx{±k + 
q/2, wi, C1J2), then using the coordinates i?, i?', r, and r', and integrating out the relative motions, we get 



By comparing this result with Eq. ([15]), we see that the 
g-integral in Eq. (jl5p is replaced with the summation 
over the Landau levels. By using the same g-cutoffs as 
in obtaining Eq. ([16]), only those states of energy levels 
el = B{2n + 1) (with = A'(0) - 2diel) in the range 
{ql/2,qll2) (in units of a = u = 1) need to be summed 
up. 

With the current-current correlation function, one then 
calculates the conductivity a according to the Kubo 
formalismi^ At very low temperatures, the correction 
to the conductivity is calculated by 

(5cr==(5P(0-,0+)/27r, (29) 

which depends on the magnetic field. The magnetocon- 
ductivity is defined as /S.(j{B) = a{B) — (t(0) where a 
is the corrected conductivity including the non-corrected 
one (see Ref. 's!) and the correction 5a. 

Shown in Fig. 5 are the present results for the magne- 
toconductivity A(T(i3) and comparison with experiments. 
For comparing with the experiments, we note that the 
overall magnitude of A(t(-B) varies largely from sample 
to sample (of the same carrier density) and from exper- 
iment to experiment i^i^ Instead to analyzing this vari- 



ation, we here confine ourselves to see the magnetic field 
dependence of Acr(-B) and therefore depict the results for 
the normalized magnetoconductivity. The solid, dashed, 
and dotted lines in Fig. 5 are the present calculations 
for the parameters {T,S,L) = (0.12 K, 8 x 10"", 1/im), 
(0.12 K, 1.5 X 10-3, l/zm), and (7 K, 2 x 10-^ 2/im), 
respectively. These sets of parameters are close to the 
conditions for the two experiments: the red filled cir- 
cles and green filled squares from Ref. [13, the up and 
down triangles and the diamonds from Ref. [2^. The 
impurity density rii = 1.15 x IQ-^a"^ and the scattering 
potential used here are the same as in the previous work^ 
that reproduces the zero-field electric conductivity of the 
experimental results?^ Clearly, the present calculation is 
in good agreement with the experimental measurements. 
Notice that there is no adjustable free parameter in the 
present calculation. 



The magnetoconductivity comes from the pseudospin 
triplet channels of the Cooperon. From Eq. ([^5]) . we see 
that the contribution SPni{B) to the conductivity from 
each Landau level of the center of mass of the Cooperon 
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r=0.12K, S= \.5xm\L=\\im 





■-^'^ . DA* 



r=7K, S=2x\0\L = 2[im 



0.0 



20 40 60 80 100 
B{mT) 



FIG. 5: (color online) Normalized magnetoconductivity as 
function of the magnetic field B. The lines are the present 
calculations: solid, dashed, and dotted lines are for the pa- 
rameter sets iT,5,L) = (0.12 K, 8 x 10"*, 1/im), (0.12 K, 
1.5 X 10~^, l/xm), and (7 K, 2 x 10""', 2/im), respectively. 
The symbols are the experimental results: red filled circles 



((5 = 8 X 10 corresponding to n = 3.15 x 10 
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green filled squares [S 
at T = 0.12 K and L 

[T = 7 K, 5 = 2 X 10 
gles (T = 4 K, (5 = 2 X 



1.5 X 10"=^ (n 
Ifim, are from Ref. 

' (n 
10" 



cm 
cm' 

1241 : up triangles 

2\ 



and 
all 



8 X 10^^ cm"^)], down trian- 
and diamonds [T = 0.26 K 



5 = 2.5 X 10"" (n = 10^^ cm"^)] all with L « 2^m are from 
Ref. H^. The magnetic field is given in unit of 10"^ Tesla. 



IV. SUMMARY 

In summary, on the basis of self-consistent Born ap- 
proximation, we have presented the WL theory of the 
Dirac fermions under the charged impurity scatterings in 
graphene. The Bethe-Salpeter matrix equations for the 
Cooperons are solved by perturbation method. There 
are three non- vanishing components in the wavefunctions 
of Cooperons under the finite-range impurity scatterings. 
This is different from the zero-range one. The pseudospin 
singlet and triplet Cooperons give rise to WL and anti- 
WL effect, respectively. For small carrier doping where 
the intervalley scatterings are much weaker than the in- 
travalley scatterings, the anti-WL effect is predominant 
over the WL one. While for large carrier doping, the WL 
effect is significant because of the intravalley scatterings 
weakened. The WL effect is also determined by the sam- 
ple size. It is found that WL is quenched at low temper- 
ature when the sample size is smaller than the inelastic 
collision length. The latter is about a few microns for 4 K 
< r < 20 K and at low carrier concentrations.^ With the 
model of charged impurity scatters justified for the elec- 
tric conductivity at zero-magnetic field, we have calcu- 
lated the magnetoconductivity. For weak external mag- 
netic field, the magnetoconductivity comes from the con- 
tributions of pseudospin triplet Cooperons. It is shown 
that the present results for the magnetoconductivity are 
in good agreement with some of the experimental mea- 
surements. 
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l-A'(0) + 2d,(2n-M)B' 



(30) 



Since A°(0) — 1, SPno is independent of B and there is 
no contribution from the pseudospin singlet channel to 
A(t{B) (both constant terms respectively in a{B) and 
(t(0) cancel each other). For the pseudospin triplet, 
l-A'(O) > 0, and at small B « [1 - A'(0)]/2d,(2n-h 1), 
5Pni{B) increases linearly with B. But the total contri- 
bution to 5a{B) does not vanish at B = because the 
Landau levels become continuum and the g-integral is re- 
stored giving rise to a constant. The total contribution 
5a{B) at _B = again cancel with the corresponding term 
in (t(0). The final result for Aa{B) depends delicately 
on B. The function Aa{B) at weak B is determined by 
the constants A'(0) and di for I ^ 0. From Fig. 3, we 
have [1- A^(0)]/4 « 2[1 - A'(0)]M|,=,,2, ee r^i. There- 
fore, the behavior of Aa{B) can be characterized by the 
constant r"^. The quantity is a measure of the in- 
tervalley scattering [because of 1 — A^ (0) as mentioned in 
Sec. III. 
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APPENDIX 

In this appendix, we discuss the impurity scattering 
potential. In our previous work^ we have illustrated how 
the intravalley and intervalley scatterings woC?) and v\ 
in Eq. ^ within the SCBA are determined from the 
microscopic electron-impurity interactions. In our nu- 
merical calculations, wo(9) and v\ are set as respectively 
the values of leading terms in their expansions for the 
charged impurities. In the derivation, the difference of 
the a and h sites in the same unit cell of graphene lattice 
was neglected. As long as the long-wave length scatter- 
ings are considered, such a difference is negligible. It is 
true for the intravalley scatterings where the long-wave 
length scatterings are predominant momentum transfers 
of the Dirac fermions. While for the intervalley scatter- 
ings, the momentum transfers are finite and more careful- 
ness are needed. Here, taking into account the sublattice 
difference, we show that only the leading terms of vq((i) 
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FIG. 6: (color online) The structure of a honeycomb lattice. 
There are two sites, a (black) and 6 (green), in each unit cell 
enclosed by the red lines. / is a vector from a to b. 



and vi need to be included in the effective potentials. 

We start with the Hamihonian for electron-impurity 
interactions in graphene, 



E 

ja 



dRua {r 



\rj-R\)n,{R) 



= ^^na{q)vai{q)n^{-q), 

qoL 



(31) 



where na{rj) is the density operator of electrons at a 
(= a or b) cite of jth unit cell of the honeycomb lattice 
(Fig. 6), ni{R) is the real space density distribution of 
impurities, and Vai{\rj — R\) is the impurity scattering 
potential. The Fourier component of the electron density 
na{q) = J2k A-qa^ka Can bc Written as 



4 

-k-qa2 



Cka2) for q^Q 



CfcQi ior q^ Q + q' ^ Q 



4-?'aiCfca2 for q 

k 



-Q 



-Q 



where k summation runs over a valley since we here 
consider only the low energy excitations. In the main 
text, the k summation ^ means and we hereafter 
use the simple notation In Eq. (j3ip . the q summation 
runs over the infinitive momentum space. By folding 



the whole space into the first Brillouin zone. Hi can be 
written as 



(32) 



kq 



where k^q run over a valley in the Brillouin zone, and 



Vdiq) Voiq - Q) 
V^{q + Q) VM 



(33) 



with 



Vd{q) 

Vo{q±Q) 

<t>a{q) = ni{-q - Qn)V{q + Qn) 

Ml) 



Ml))' Mq) 

(l)a{q±Q)atc 
(f)b{q±Q)atc 



ni{-q - Qn)V{q + Qn)- 



Here V{q) is the Fourier component of Vai{r)^ and Qn 
is the reciprocal lattice vector. There is a phase factor 
g»(9+Qi.) ' the expansion of 4>b{q) because of the posi- 
tion difference I between a and b sites. The expressions 
for 4>a,b[q)atc are similar to 4>a,b{q) except avoiding triple 
counting since 4>a,biq ± Q)atc imply the intervalley scat- 
terings. [There are other two equivalent valleys for each 
valley indicated in Fig. 1(a).] For example, the leading 
oi-der in ^q - Q)atc ~ 'l>a{-Q)atc is ni{Q)V{Q) where 
Q (with IQI = 47r/3a, a ~ 2.4 A as the lattice constant) 
is a momentum difference between the nearest-neighbor 
Dirac points in the Brillouin zone. For this Q, there are 
other two vectors of the same magnitude differ from Q by 
two reciprocal lattice vectors [see Fig. 1(a)], respectively. 
These two terms should be excluded from the summation. 
Since the case of q is under consideration, we approx- 
imate the off-diagonal potential as Vo{q ±Q) ~ Vo{±Q)- 
We have argued^ that if the effective impurity scatter- 
ings are given as Eq. then the vertex Tx{k,uji,uj2) 
can be expanded as Eq. ([5])- [The off-diagonal parts 
in Eq. ^ are different from that in Ref. ^ where 
ctq was used instead of ci because the basis was Tp], = 
(Cfcai' 4a2' 4b2) '^^^^ reflected y-axis in valley 2. 
But the result for Tx{k,uJi,uj2) is unchanged.] Here, sup- 
pose the vertex Tx{k, wi, W2) is given as Eq. ([5]), we want 
to see how ^0(9) and vi in Eq. ([2] are expected from 
the microscopic potentials given by Eq. ([33]). The x- 
direction current vertex vTx{k^uJi^uj2) under the SCBA 
satisfies the following 4x4 matrix integral equation, 
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T,{k,LUi, UJ2) = T3a, + ^J2^V^{k-k')G{k',0Ji)T,{k\uJi,LU2)G{k\uj2)V,{k' - fc)), (34) 

k' 

where (• • • ) means the average over the hnpurity distributions [Fig. 1(c)]. Notice that the product 
G{k,u}i)rx{k,uji,u}2)G{k,uj2) can be expanded in A^{k), 

G{k,uji)Tx{k,uJi,uj2)G{k,uj2) = A^{k)Ljj> (k, uji,uj2)vf (fc, 0.^1, 0.^2), 

3 3' 

where the functions Ljj' {k,uji,uj2) have been defined in Ref. S We then need to calculate the expectations {Vi{k — 
k')A^{k')V,{k' - k)) in Eq. m- Firstly, we calculate {V,{q)A^V,{-q)) = A^[{M'l)M-<l)) " {MQ)atcM-Q)atc)- 
Using the expressions for <pa,biq), we obtain 

n 

where the use of the fact that J2n^o l^il + QrOpe^*'^" ' = for g ^ has been made. Similarly, we get 
{MQ)atcM-Q)atc) = Vn,Y,'VHQ + Qn)e-'(^+'^-'>-'^- Vn.V^Q), 

n 

where J2n ' means avoiding triple counting, and Q ■ ^ = because for I there is a Q orthogonal to it [see Figs. 1(a) 
and 5]. For Aj{k') with j ^ 0, we note 

^^ 2(^0 = A^.2(fc)(cos6'±icr^sin6l), Al[k!) = (fc) (cos 26* + icr^ sin 26*), (35) 

where B is the angle between k and k' . Since V[k' — fc) depends on B through cos 6*, the integrals of the integrands 
with factor sin 6* or sin 26* vanish. (V^i(fc — k')A^ {k')Vi{k' — k)) for j 7^ arc then calculated as 

{Vi{k'^k')Al2ik')Vi{k' - k)) -> Al^Vni^V^Hk ~ k' + Q„|)cos6i w Al^VuiV^Hk/ - k\)coaB 

n 

{Vi{k - k')A^{k')V^{k' ~ k)) -> A^Vn,"Yv^{\k - k' + Qn\)e-'^'^+'^'^^-''cos29 = A^Vn.,V^{\k' ~ fc|)cos26'. 

n 

By defining VQ{q) = V'^{q) and — V'^{Q), we obtain exactly the same equation as Eq. (11) in Ref. 3 for determining 
yj(k,uji,oj2)- Thus, we have proved that only the leading terms of VQ{q) and vi are necessarily taken into account in 
the current vertex corrections. 

On the other hand, under the SCBA, the self-energy is given by 

S(fc,07) = :^^(V^,(fc-fc')G(fc',C^)V^^(fc'-fc)) 

k' 

« y E^E ^'(^ - ^' + Q") + - g„)]Go(fc',c^) + V\k - k')G,{k' ,uj)k ■ k'k ■ arj 

A;' n n 

~ y E{["^'(^ " + y^(Q)]Go{k',L^) + V^ik^ k')G,{k\uj)k ■ k'k ■ or,}, 

k' 

I 

which is consistent with Eqs. ([3|) and ([4]). Therefore, the Eq. ([2|) with wo(q) and v\ defined above, 
theory is simplified with the effective potential given by 
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